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The Black Ring is Unstable 
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We study non-axisymmetric linearised gravitational perturbations of the Emparan-Reall black ring 
using numerical methods. We find an unstable mode whose onset lies within the “fat” branch of 
the black ring and continues into the “thin” branch. Together with previous results using Penrose 
inequalities that fat black rings are unstable, this provides numerical evidence that the entire black 
ring family is unstable. 


Introduction. During the golden age of general rel¬ 
ativity, four-dimensional black holes were shown to be 
remarkably featureless. The Kerr black hole is the most 
general asymptotically flat black hole solution to the vac¬ 
uum Einstein equation pQ . It has spherical topology and 
depends only on two parameters: the mass M and the 
angular momentum J mm- Moreover, for specific val¬ 
ues of M and J, there is only one such black hole, i.e. 
the Kerr black hole is unique. Additionally, its linear 
stability has been shown for all modes [3]. Numerical 
simulations suggest that Kerr black holes are also non- 
linearly stable, though a complete proof remains illusive 
(see however [5HZ]). These facts have lead to the conjec¬ 
ture that ‘four-dimensional black holes have no hair’. 

There is a natural higher-dimensional extension of the 
Kerr black hole - the Myers-Perry black hole [5]. Like 
Kerr, it has spherical topology, depends uniquely on the 
mass and angular momenta, and (for sufficiently slow 
rotation) has substantial evidence for its stability [9, |i TD], 

Nevertheless, after the initial discovery of the black 
ring by Emparan and Reall mm, it became clear that 
black holes in higher dimensions are dramatically dif¬ 
ferent from those in four-dimensions. These black rings 
are five-dimensional solutions with topology S 1 x S 2 and 
come in two types described by their shape: “fat” rings 
and “thin” rings. The two branches of solutions are dis¬ 
tinct except for a unique ring that is both fat and thin. 
For every fat ring, there is a thin ring and a Myers-Perry 
black hole with the same mass and angular momentum. 
Thus, uniqueness is broken among black rings, but also 
with Myers-Perry black holes. In light of this, a higher¬ 
dimensional no-hair theorem seems unlikely to hold. 

Yet, another surprise in higher dimensions is the exis¬ 
tence of a new type of instability. The so-called Gregory- 
Laflamme instability was first found in black strings |13| , 
but there are analogous instabilities for rapidly rotating 
Myers-Perry black holes mmm- The general picture 
is that black objects with extended directions are unsta¬ 
ble to perturbations along those extended directions. 

Since many of the solutions that violate uniqueness 
have extended directions, there is hope that the spirit of 
the no-hair theorems can be restored. That is, there is a 
unique stable solution which is a slowly rotating Myers- 
Perry black hole - a dynamical no-hair conjecture. 

In five dimensions, the (in)stability of the black ring 


is a natural setting to test this conjecture. Indeed, they 
violate uniqueness, but moreover all known solutions in 
five dimensions are either Myers-Perry or contain a topo¬ 
logically S 1 x S 2 black object as a horizon component. 

Fat rings were conjectured to be unstable in mm 
and confirmed in l20j using local Penrose inequalities. 
However, little is known about the stability of thin rings, 
except that very thin rings ought to be unstable to 
Gregory-Laflamme because they resemble boosted black 
strings EHEa EH- Though this argument is physi¬ 
cally sound, it says little about the stability of thin, 
but not very thin, rings (which happens to coincide with 
where black rings violate uniqueness with themselves and 
Myers-Perry.) Our goal is to demonstrate that such a 
window of stability does not exist by computing gravita¬ 
tional perturbations of the black ring. 

Numerical Approach. Before embarking on a study 
of gravitational perturbations, let us present coordinates 
for the black ring that we found suitable for numerics: 


d s 2 =R Z l — (1 — y 2 ) 2 F dt 


V 

h 4 


4 ^ dy 2 + y 2 (2~y 2 )S(d'ip-nh 4 Wdt) 2 


2-y 2 g 


4 ^ 2 f 2 , o 2„,2 


2 - z 2 h 


dx 2 + ff z x 2 {2 — ar)(l — x 2 ) 2 fi d <f> 


(1) 


where the functions /i, f 2 , g , h, F, S , and W are 
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and a > 1, (3 > 0 are constants that parametrise this 
family of solutions (the constant R merely sets a scale). 
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These parameters are related to the more familiar ones 
presented in [12] via 


/3 2 x a 2 (3 2 
2 + /3 2 ’ ~ 2 + a 2 /? 2 


(4) 


which satisfy 0 < v < A < 1. The parameters A and v 
are more prevalent in the literature, so we will refer to 
these rather than a and ft. The remainder of the line 
element in [121 can be reproduced by the redefinitions 
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where quantities with a tilde refer to those in |12| . 

Our coordinates range in x £ [0,1] and y £ [0,1] with 
the axis of rotation at y = 0, the horizon at y = 1, 
the outer axis of the ring at x = 0, and the inner axis 
at x = 1 . Asymptotic infinity is at x = y = 0, which 
corresponds to h = 0. The period of ^ and (j) are set to 
27 t and the temperature of the horizon is 1/ (2n). 

Unless A = 2^/(1 + v 2 ), there will be a conical singu¬ 
larity at the inner axis. We shall see that the singular 
solutions will be useful to us. In particular, there is a 
static solution when A = v. We are chiefly interested 
in the non-singular family of solutions which we call the 
“balanced” ring. When the ring is balanced, 0 < v < 1/2 
are the thin rings, while 1/2 < v < 1 are the fat rings. 

We wish to study gravitational perturbations of ([!]). 
Let us write the perturbed metric as g a b = °5ab + h a b , 
where °g a b refers to the background solution 0. and hab 
is our metric perturbation. Since we are working with 
the vacuum Einstein equation, we are free to impose the 
transverse-traceless gauge condition 


X7 a h ab = 0 , h a a = 0 , (6) 


where, as elsewhere in this manuscript, covariant differ¬ 
entiation and the raising and lowering of indices are done 
with respect to the background °g a b- In this gauge, the 
linearised Einstein equation is 

(A L h)ab = -V c V c /i ab - 2 R a c b d h cd = 0 , (7) 


where is the Lichnerowicz operator. 

Thin rings are extended in the d^, direction, so the 
Gregory-Laflamme instability in the black ring would 
break the d^ symmetry. For simplicity, we preserve the 
remaining symmetry d We therefore perform a mode 
decomposition h a b = e~ lu)t+lm '‘ lj h a b , where h a b are func¬ 
tions of x and y. 

Our background has a fixed temperature T = l/(27r), 
so the frequency w is equivalent to the more general di¬ 
mensionless quantity vj = lo/(2ttT). Incidentally, this 


quantity is equivalent to w = (i? 0 ut — i?i n )w, where R ou t 
and i?i„ are the outer and inner S 1 equatorial horizon 
radii, respectively. 

Since the to = 0 modes do not break the rotation axis 
and to = 1 modes were not found to be unstable in other 
systems with non-axisymmetric instabilities m, we will 
for definitiveness and simplicity set m = 2. 

Preserving d$ symmetry lets us set h^ = 0 for /z / <j>. 
This leaves 11 functions, one of which can be removed 
by imposing tracelessness. After imposing tracelessness, 
there are six components of 0 which together with the 
four non-trivial components of the transverse condition 
in 0 form a set of 10 independent equations. The re¬ 
maining components of ([T]) can be derived from this set. 

Our task is to solve this set of 10, two-dimensional 
partial differential equations in the form of a quadratic 
eigenvalue problem in ui. As boundary conditions, we 
impose regularity on the outer axis, ingoing boundary 
conditions at the horizon, and outgoing boundary condi¬ 
tions at infinity. On the inner axis, we demand that the 
conical excess/deficit does not change, which is equiva¬ 
lent to regularity when there is no conical singularity. 

This problem is complicated by the fact that infinity 
is at the coordinate singularity x = y = 0. Because 
of this singularity, approaching the point x = y = 0 
from different directions will yield different values for h a b- 
To obtain a well-posed problem, our remedy is to use a 
different coordinate system near infinity given by 


7 . r px^/2—3^ 

P = h , £ = y 1-^-, (8) 

where the function h was given in j2j|. In these coor¬ 
dinates, spatial infinity is at the hyperslice p = 0. We 
divide our domain into two non-overlapping coordinate 
patches, one in (p, £) coordinates containing infinity, and 
another in (x, y) coordinates containing the horizon and 
the inner axis. Grids can then be placed on these patches 
using transfinite interpolation. We must also impose ad¬ 
ditional patching conditions that require that h a b and its 
first derivatives match on patch boundaries. 

Once the problem is discretised by such a patched grid 
(we use pseudospectral collocation on Chebsysliev grids), 
it is reduced to a quadratic eigenvalue problem in linear 
algebra, which can be solved on a computer. Unfortu¬ 
nately, because of the size of the matrices required, a di¬ 
rect computation of the spectrum (reduction to a linear 
matrix pencil followed by QZ factorisation) yields a large 
number of non-physical spurious modes. We were unable 
to extract the physical modes from such a spectrum. 

We must therefore seek an alternate route. Rather 
than the entire spectrum, we wish only to obtain the 
unstable Gregory-Laflamme mode of the balanced ring. 
By varying the parameters of the black ring, it is natural 
to suspect that this mode is connected to some Gregory- 
Laflamme mode of the static ring, the onset of which is 
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a zero frequency mode (ui = 0). 

Since we suspect that our desired mode is connected to 
a zero frequency mode, consider the more general prob¬ 
lem (with the same gauge and boundary conditions) 

A L h ab = -k 2 h ab . (9) 

This problem arises in perturbations of the six- 
dimensional solution (black ring) x R. If there is a zero 
frequency mode solution to A Lh a b = 0, then there is 
also a zero frequency mode solution to © with k = 0. 
One can in principle find this mode by setting lo = 0, 
and solving © for k while varying the parameters until 
k = 0. Indeed, this method has been used successfully in 
the past puns]. 

For our purposes, we need not find the zero frequency 
mode of the static ring. We merely observe that there 
are solutions of © on the static ring with uj = 0 and 
k ^ 0 that are connected to the zero frequency mode 
of the static ring, and are hence also connected to our 
desired Gregory-Laflamme mode of the balanced ring. 

Given this, we first solve ([9]) with w = 0 on the static 
ring for k. This problem is substantially simpler than 
solving |7| directly. This is a linear eigenvalue problem 
in k 2 . The static background and ui = 0 introduces ex¬ 
tra symmetry and reduces the number of functions from 
10 to 7. A suitable definition of the perturbation func¬ 
tions also yields matrices for the eigenvalue problem that 
are purely real. Furthermore, we know (as we know for 
black strings) that the real, positive k 2 modes are con¬ 
nected to the zero frequency mode, significantly reduc¬ 
ing our search space for physical modes. Attempting this 
problem for a particular static ring by QZ factorisation 
with several grid resolutions yielded a single real, positive 
value of k 2 . 

Having obtained a single solution connected to the so¬ 
lutions we are after, we can proceed by repeated appli¬ 
cation of Newton-Raphson on several grid resolutions. 
Since the Gregory-Laflamme modes of the static ring are 
pure imaginary, we can increase T = —iui until we find 
that k = 0. This puts us on a solution of ([7|. From here, 
we can solve for oj while varying the ring parameters. 
We increase the rotation until the ring is balanced and 
we have thus arrived on a desired solution. We hence¬ 
forth take the black ring to be balanced and present our 
results as we vary the parameter v. 

As a check, we have explicitly evaluated the Geroch- 
Held-Penrose scalars constructed in [22] and found non¬ 
zero values, confirming that this mode cannot be pure 
gauge. More details of our calculation and a number of 
numerical checks can be found in the appendix. 

Results. FIG.[l]and[2]present the imaginary and real 
parts, respectively, of the dimensionless quantity vo = 
uj/(2nT) as a function of v. (Recall that thin rings have 
0 < v < 1/2 and fat rings have 1/2 < v < 1.) Dots 
and squares represent different resolutions, where ( N + 
N) x N refers to N 2 points per patch. The agreement 



FIG. 1: Imaginary part of the frequency vj = oj/(2tyT) for 
non-axisymmetric perturbations with m = 2 as a function of 
v at two resolutions. The coloured region corresponds to fat 
rings which have previously been shown to be unstable. The 
solid line is a polynomial extrapolation. 



FIG. 2: Real part of the frequency vj = co/(2nT) as a function 
of v. Same colour scheme as FIG. [T] 

between resolutions is reassuring. Data points below v ~ 
0.144 and above v ~ 0.52 were discarded due to lack of 
numerical precision. (Our criterion is to discard values 
where the two resolutions differ by more than 0.1%). The 
vertical dashed line indicates the separation between fat 
and thin rings, with the coloured region corresponding to 
fat rings. 
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FIG. 3: Isometric embedding of constant ip slices of the S 2 
spatial sections of the perturbed black ring horizon (see ap¬ 
pendix for more details). This plot corresponds to v = 0.2. 

Our main result can be seen in Fig. [l] where positive 
Im(G7) indicates an instability. An instability seems to 
exist for all values of v < 0.55 (the upper bound is an 
extrapolation of data with v > 0.3 using a second order 
polynomial in is). In particular, this instability extends 
into a region of the fat rings with 0.5 < is < 0.55. Since 
all fat rings are already unstable to axisymmetric pertur¬ 
bations rrSl - DHTl , this would imply that the Emparan-Reall 
ring is unstable for all ranges of parameters. We stress 
that though we do not have points below v < 0.144, the 
Gregory-Laflamme argument should be valid for small is, 
so we expect the instability to persist down to is = 0. 

We note that the local Penrose inequalities suggest 
that the axisymmetric (m = 0) instability should be 
marginal at v = 1/2, so fat rings near v = 1/2 should 
have a small growth rate. Our results therefore suggest 
that the in = 2 mode is dominant over the m = 0 inodes 
for fat rings near is = 1/2. For larger values of is, the 
m = 2 modes are no longer unstable, so the m = 0 modes 
should dominate. Finding the transition point would re¬ 
quires data on the m = 0 sector of perturbations, which 
we leave for future work. 

A curious fact about FIG. [2] involves the so-called su¬ 
perradiant bound. Based on the results of [ 221 12 H , the 
onset of the instability must satisfy 0 < Re(w) < mf2, as 
we have verified. Away from the onset at v < 0.305, we 
find instead that Re(w) > mil. Though this is not in con¬ 
flict with |23l I24| , we do not have many examples where 
Im(w) > 0 and Re(w) > mfl. The physical significance 
of this property is unclear at the moment. 

To visualise the horizon evolution under this instabil¬ 
ity, we have constructed an isometric embedding for our 
perturbations (see appendix for more details). The re¬ 


sults are depicted in FIG. |3J and are similar to those ob¬ 
tained in [25] for the Gregory-Laflamme instability. This 
corroborates our claim that these non-axisymmetric per¬ 
turbations are similar to that of the black string. The 
general shape of the curve in Fig. [I] is also characteristic 
of Gregory-Laflamme inodes mm- 

Outlook. We have given substantial numerical evi¬ 
dence that the entire Emparan-Reall black ring is un¬ 
stable. Fat rings are unstable to axisymmetric pertur¬ 
bations, and thin rings are unstable to non-axisymmetric 
perturbations resembling the Gregory-Laflamme instabil¬ 
ity of the black string. There is a competition between 
these instabilities for fat rings near is = 1/2. 

We focused our computation entirely on the m = 2 
inodes. It would be interesting to see how the other 
modes behave. A study of the m = 0 inodes would eluci¬ 
date its competition with the to = 2 modes. There could 
also be instabilities corresponding to to = 1 , though 
these were absent in Myers-Perry [TO]. Perturbations 
with to > 2 correspond to shorter wavelength and may 
have an onset for smaller is than that of to = 2. One can 
also consider perturbations that break the d$ symmetry. 

While we have only studied stability of the Emparan- 
Reall black ring, there are many other solutions. In five 
dimensions, the black ring can also rotate in the d$ di¬ 
rection as in the double-spinning ring [26] and the helical 
rings m- A study of rotating black strings [T5] suggests 
that the instability in doubly-spinning rings would have 
a higher growth rate. There are also multi-horizon solu¬ 
tions ( e.g. black Saturns [28], di-rings (2j.il■ and bicycling 
rings BDj). These solutions contain black rings as horizon 
components and might share many of the same stability 
properties. In six and higher dimensions, there are also 
black rings [SI Eli and their associated multi-horizon so¬ 
lutions, but there are additionally ringoids [271133] and 
lumpy black holes [32] El] ■ Little is known about the sta¬ 
bility of these solutions, though many of them resemble 
an unstable Myers-Perry black hole, or a black ring, and 
hence might also share many stability properties. The 
addition of matter may add a stabilising effect, particu¬ 
larly in supersymmetric setups [35ti38j . 

The endpoint of these instabilities remains an impor¬ 
tant and open problem. Work is in progress [32] ■ Ax¬ 
isymmetric instabilities are expected to lead towards 
Myers-Perry black holes. For non-axisymmetric pertur¬ 
bations of very thin rings, the Gregory-Laflamme insta¬ 
bility in the (unboosted) black string suggests that black 
rings would develop a naked singularity and violate cos¬ 
mic censorship [40]. The similarity of our results to the 
Gregory-Laflamme instability seems to support this idea. 
It is important to note, however, that this instability 
in the black ring emits gravitational radiation. Even 
for very thin rings, the solutions resemble boosted black 
strings, and it is unclear what role this would play in the 
overall time evolution. 
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Appendix 

Numerical Details. Here, we supply more details to our numerical approach that were absent in the main text. 
Our coordinates for the background 0 are not the usual coordinates used to describe black rings such as that of 
m It was chosen for numerical convenience (for example, many of of the boundary conditions are simpler in these 
coordinates). In these coordinates, the expressions for the mass (M), angular momentum (J), temperature (T), 
horizon angular frequency (H), and horizon area (A) are 


M = 

J = 

T = 

n = 

A = 
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( 10 ) 


where G is the gravitational constant and we have temporarily rescaled t —» t/R to more clearly show how these 
quantities depend on the dimensionful scale R. The quantities to the right of the arrows refer to the balanced ring. 
Now we write down an ansatz for our perturbed metric in full: 
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where the SH J s are unknown functions of x and y , and the remaining functions and constants have been defined 
previously in (|2j) and ([3]). This metric is already of the form g a b = °g a b + h a b, where °g a b is the background solution 
and the SH’s appear linearly in h a b- Now we remove the factors 


. _ g — iut+imit e ik 0 u/h^3/2^ _ y2^—i(w—mQ.) 


yV 2 — 


SH, , 


( 12 ) 


and work with the SH ’s instead. The factor of j S) 0 f course, our mode decomposition of the perturbations. 

The remaining factors and some of the extra factors in (11) correspond to the behaviour of the perturbation functions 
after imposing the appropriate boundary conditions. These factors can be derived using an approach similar to that 
of P]3]. Note that some of these factors depend on m > 2 and will change for m = 0 and m = 1. 
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The equations we wish to solve are 


h U a =0 


y a hab = 0 , A Lhab = 0 . 


(13) 


We can use the traceless equation to eliminate 6H 5 from the remaining equations. This leaves us with 10 functions to 
solve for. The linearised Einstein equation A Lh ab = 0 produces second-order equations in the SH’s of which we only 
keep those associated with SHj for j £ {2,3, 7, 8, 9,11}. Together with the transverse condition V“/i 0 & = 0, this forms 
our set of 10 equations that we will solve numerically. One can verify that the remaining components of A Lh ab = 0 
can be derived from this set of equations and its derivatives. Because of the coordinate singularity at x = y = 0, we 
must also obtain this same set of equations in the (p, £) coordinate system defined in <©• 

Now we discuss boundary conditions. The factors we removed in (12) and our choice of ansatz (11) guarantee 


that our boundary conditions are satisfied so long as the SH's remain finite. Nevertheless, for numerical accuracy we 
impose the following boundary conditions, all of which can be derived from a series expansion of our equations. At 
x = 0, we impose 


SHi + SH 2 + (1 - y 2 f{25H 3 + 5H A ) = 0 , 
d x 6H 3 = 0 , j € (1,2,4,6,7,8,9,10,11} . 

At x = 1, we impose 

y 2 (2 - Z/ 2 )Affi + (y 2 (2 - y 2 ) + P 2 )SH 2 + (1 - y 2 ) 2 (2y 2 (2 - y 2 )5H 3 + (y 2 (2 - y 2 ) + /3 2 )8H 4 ) = 0 , 

d x 5H 3 =0, j e {1,2,4,6,7,8,9,10,11} . 


At y = 0, we impose 


5H 2 + SH 4 = 0 , 

SH 7 - 26Hu = 0 , 

6H S - 26H 6 = 0 , 

SH 10 + 2f 2 6H 4 = 0 , 
dySHj = 0 , j e {1,3,4,6,9,11} . 


And at y = 1, we impose 


5H X + f\5H 2 = 0 , 


6H 2 - 
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SHs = 0 


6H 7 - 2 v / 1 + x 2 (2-x 2 )/3 2 5H 9 = 0 , 
6H 10 - 2i/l + £ 2 (2 — x 2 )/3 2 6H 6 = 0 , 


dySHj = 0 , j G {3,4,6,8,9,11} . 


(14) 


(15) 


(16) 


(17) 


Under the coordinate transformation ([8]) there are boundary conditions at £ = 0 and £ = 1 that are equivalent to those 
at y = 0 and x = 0, respectively. At infinity (p = 0), we impose the lowest-order series expansion of the equations 
about p = 0. There, the boundary conditions are too long to merit their inclusion in this manuscript. 

As mentioned in the main text, we solve the equations above using Newton-Raphson. As a means of obtaining a 
seed, we set the background to the static ring ( i.e. a = 1) and solve 


k 2 

h a a = 0, V a hab = 0, A L h ab = —k 2 h ab = - 2 h ab . (18) 

i\ z Kq 

We can take the same set of equations and boundary conditions we had before with one difference. When k 0, the 
functions near infinity (p = 0) go as SH ~ e~ k ^ p , so there we impose the Dirichlet conditions SHj = 0. 

Since we are on the static ring, we have fl = 0. If we further redefine u = iT, one can show that all factors of i 
drop out of the equations of motion and boundary conditions. Hence, we are solving an eigenvalue problem with real 
matrices. 
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To obtain our first solution, we set T = 0 and view (18) as a linear eigenvalue problem in k 2 . (Note that had we 


factored out e from the 8 H' s, we would have been left instead with a quadratic eigenvalue problem in k.) In this 
case, there are extra symmetries that allow us to set SHq = 8 H$ = 5Hg = 0. For a particular static ring (/3 = 1 was 
a good choice), we were able to find a single positive, real value of k 2 via QZ factorisation. 

Having obtained our first solution, we continue to view (18) as a linear eigenvalue problem in k 2 , but now solve 


it using Newton-Raphson. We slowly increase T until we find that k is close to zero. Note that if T ^ 0, we are no 
longer allowed to set SH e = 8 H$ = 5Hg = 0, and must solve for the full set of equations. 


Once k is close to zero, we can input the result as a seed for (13), viewed as a quadratic eigenvalue problem in F. 


We then increase a (and hence also f l) until the ring is balanced. Note that when fi / 0, we expect T to now be 
complex. Once balanced, we can explore the rest of the parameter space of the balanced ring. 

At this point, we give a few numerical results to our calculation. A plot of the eigenfunction 8 H 4 is given in FIG. [4] 
The apparent smoothness of the function is reassuring. From these plots, one can see the non-overlapping patches we 
have placed with transfinite interpolation. One can also see the apparent coordinate singularity at x = y = 0 and that 
it is removed when mapped to p = 0 in the (p, £) coordinate system. We have chosen to normalise the eigenfunctions 
so that SH 1 = 1 + i at the corner where the horizon meets the inner axis. 

We have performed a convergence test, the results of which are displayed in Fig. [5j We monitored the quantity 
|1 — w/v/ wn+i| for several different resolutions. Here, lon is the complex frequency of the unstable mode with m = 2 
described in the main text, computed for v = 1/2. The results are plotted in a log scale in Fig. [ 5 J where a straight 
line can be seen. This is good evidence that the convergence of our method is exponential, as dictated by spectral 
collocation methods. 

Embedding Diagrams. Now let us construct an embedding diagram. Recall that the ring has horizon topology 
S 1 x S 2 . We will embed the S 2 of this horizon into R 3 , and plot the result as a function of the S 1 coordinate. This 
partially follows [Hi- 

First, let us move to ingoing Eddington-Finkelstein coordinates 


= t + log(l - y 2 ) , $ = ip + n\og(l - y 2 ) . 


(19) 


In these coordinates, we have 




(i -y 1 ) 


2_ ^-iujv-\-imip 


( 20 ) 


so constant v and ip slices on the horizon have nontrivial perturbations. At constant v and 1 p, the induced metric on 
the horizon takes the form 


gxxixP)dx T Q(fxj){sx)d(p . (21) 

Our aim is to embed the S 2 of the ring horizon into R 3 : 

dX 2 + dY 2 + Y 2 d(p 2 , (22) 

which we have chosen to write in cylindrical coordinates. Now let us suppose we have a parametrised curve X(x) and 
Y(x). Then the induced metric on this curve is 

(. X’(x ) 2 + Y'{xf)dx 2 + Y {x) 2 (p 2 . (23) 


Equating this with (21) gives 


Y ( x ) = \Jgu{x) , X'{x) = \/g xx {x) - Y'(x ) 2 , 


(24) 


which is an ordinary differential equation we can solve to compute the embedding diagram. 

To fix units and normalisations, we rescale GM = 1, fix the Eddington-Finkelstein coordinate v = 0, and normalise 
the eigenfunctions so that 8 H$ = 1 where the horizon meets the inner axis. The integration constant above is fixed so 
that the S 2 is centred on the coordinates ( i.e. the maximum and minimum value of X on the embedding differ only 
by a sign), with the inner axis of the horizon at negative A' and the outer axis at positive X. Though this embedding 
is in R 3 , the symmetry about <p lets us ignore this coordinate an plot just X(x) and Y(x). We do so as a function of 
-i/c The results of this embedding are shown in FIG. [6] for two values of v. 









FIG. 4: Real (top plots) and Imaginary (bottom plots) parts of 8 H 4 in the (x, y ) coordinate system (left plots) and (p, £) 
coordinate system (right plots). 
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FIG. 5: Convergence test as a function of the number of points N. The resolutions used were of the form (N + N) x N, 
indicating N 2 points per patch. 




FIG. 6: Isometric embedding of constant ip slices of the S 2 spatial sections of the perturbed black ring horizon. The left plot 
corresponds to v = 0.5, while the right corresponds to v = 0.2. 


* Electronic address: J.E.Santos@damtp.cam.ac.uk 
t Electronic address: B.Way@damtp.cam.ac.uk 

[1] R. P. Kerr, Gravitational field of a spinning mass as an 
example of algebraically special metrics, Phys. Rev. Lett. 
11 (Sep, 1963) 237-238. 

[2] D. Robinson, Four decades of black holes uniqueness 
theorems, in The Kerr Spacetime: Rotating Black Holes 
in General Relativity (D. L. Wiltshire, M. Visser, S. M. 
Scott, ed.). Cambridge University Press, 2009. 

[3] P. T. Chrusciel, J. L. Costa, and M. Heusler, Stationary 
Black Holes: Uniqueness and Beyond, Living Rev.Rel. 


15 (2012) 7, arXiv: 1205.6112. 

[4] B. F. Whiting, Mode Stability of the Kerr Black Hole, 

J.Math.Phys. 30 (1989) 1301. 

[5] M. Dafermos and I. Rodnianski, Decay for solutions of 
the wave equation on Kerr exterior spacetimes I-II: The 
cases |a| <C M or axisymmetry, arXiv: 1010.5132 

[6] M. Dafermos, I. Rodnianski, and 

Y. Shlapentokh-Rothman, Decay for solutions of the 
wave equation on Kerr exterior spacetimes III: The full 
subextremal case |a| < M, arXiv: 1402.7034 

[7] M. Dafermos, I. Rodnianski, and 

Y. Shlapentokh-Rothman, A scattering theory for the 
wave equation on Kerr black hole exteriors, 
arXiv:1412.8379 







10 


[8] R. C. Myers and M. Perry, Black Holes in Higher 
Dimensional Space-Times, Annals Phys. 172 (1986) 

304. 

[9] M. Shibata and H. Yoshino, Bar-mode instability of 
rapidly spinning black hole in higher dimensions: 
Numerical simulation in general relativity, Phys. Rev. 
D81 (2010) 104035, arXiv: 1004.4970 . 

[10] O. J. Dias, G. S. Hartnett, and J. E. Santos, 
Quasinormal modes of asymptotically flat rotating black 
holes, Class.Quant.Grav. 31 (2014), no. 24 245011, 
arXiv:1402.7047 . 

[11] R. Emparan and H. S. Reall, A Rotating black ring 
solution in five-dimensions, Phys. Rev. Lett. 88 (2002) 
101101, hep-th/0110260 . 

[12] R. Emparan and H. S. Reall, Black Rings, 

Class. Quant. Grav. 23 (2006) R169, hep-th/0608012 . 

[13] R. Gregory and R. Laflamme, Black strings and 
p-branes are unstable, Phys.Rev.Lett. 70 (1993) 
2837-2840, hep-th/9301052 . 

[14] O. J. Dias, P. Figueras, R. Monteiro, J. E. Santos, and 
R. Emparan, Instability and new phases of 
higher-dimensional rotating black holes, Phys.Rev. D80 
(2009) 111701, arXiv:0907.2248 . 

[15] O. J. Dias, P. Figueras, R. Monteiro, H. S. Reall, and 
J. E. Santos, An instability of higher-dimensional 
rotating black holes, JHEP 1005 (2010) 076, 
arXiv:1001.4527. 

[16] M. Shibata and H. Yoshino, Nonaxisymmetric instability 
of rapidly rotating black hole in five dimensions, 

Phys.Rev. D81 (2010) 021501, arXiv:0912.3606. 

[17] G. S. Hartnett and J. E. Santos, Non-Axisymmetric 
Instability of Rotating Black Holes in Higher 
Dimensions, Phys.Rev. D88 (2013) 041505, 
arXiv:1306.4318 . 

[18] G. Arcioni and E. Lozano-Tellechea, Stability and 
critical phenomena of black holes and black rings, 

Phys.Rev. D72 (2005) 104021, hep-th/0412118. 

[19] H. Elvang, R. Emparan, and A. Virmani, Dynamics and 
stability of black rings, JHEP 0612 (2006) 074, 
hep-th/0608076 . 

[20] P. Figueras, K. Murata, and H. S. Reall, Black hole 
instabilities and local Penrose inequalities, 

Class. Quant. Grav. 28 (2011) 225030, 

arXiv:1107.5785 . 

[21] J. Hovdebo and R. C. Myers, Black rings, boosted 
strings and Gregory-Laflamme, Phys.Rev. D73 (2006) 
084013, hep-th/0601079 . 

[22] M. Godazgar and H. S. Reall, Peeling of the Weyl 
tensor and gravitational radiation in higher dimensions, 
Phys.Rev. D85 (2012) 084021, arXiv: 1201.4373 . 

[23] S. Teukolsky and W. Press, Perturbations of a rotating 


black hole. Ill - Interaction of the hole with 
gravitational and electromagnet ic radiation, 

Astrophys.J. 193 (1974) 443-461. 

[24] R. M. Wald, General Relativity. Chicago University 
Press, Chicago, USA, 1984. Chapter 12. 

[25] R. Gregory and R. Laflamme, The Instability of charged 
black strings and p-branes, Nucl.Phys. B428 (1994) 
399-434, hep-th/9404071 . 

[26] A. Pomeransky and R. Sen’kov, Black ring with two 
angular momenta, hep-th/0612005 

[27] R. Emparan, T. Harmark, V. Niarchos, and N. A. 

Obers, New Horizons for Black Holes and Branes , 

JHEP 1004 (2010) 046, arXiv:0912.2352 . 

[28] H. Elvang and P. Figueras, Black Saturn, JHEP 0705 
(2007) 050, hep-th/0701035 . 

[29] H. Iguchi and T. Mishima, Black di-ring and infinite 
nonuniqueness, Phys.Rev. D75 (2007) 064018, 
hep-th/0701043 . 

[30] H. Elvang and M. J. Rodriguez, Bicycling Black Rings, 
JHEP 0804 (2008) 045, arXiv:0712.2425 . 

[31] B. Kleihaus, J. Kunz, and E. Radu, Black rings in six 
dimensions, Phys.Lett. B718 (2013) 1073-1077, 

arXiv:1205.5437 . 

[32] O. J. Dias, J. E. Santos, and B. Way, Rings, Ripples, 
and Rotation: Connecting Black Holes to Black Rings, 
JHEP 1407 (2014) 045, arXiv: 1402.6345 . 

[33] B. Kleihaus, J. Kunz, and E. Radu, Black ringoids: 
spinning balanced black objects in d> 5 dimensions - 
the codimension-two case, JHEP 1501 (2015) 117, 
arXiv:1410.0581. 

[34] R. Emparan, P. Figueras, and M. Martinez, Bumpy 
black holes, JHEP 1412 (2014) 072, arXiv: 1410.4764 . 

[35] H. Elvang, R. Emparan, D. Mateos, and H. S. Reall, A 
Supersymmetric black ring, Phys.Rev.Lett. 93 (2004) 
211302, hep-th/0407065 . 

[36] I. Bena and N. P. Warner, One ring to rule them all ... 
and in the darkness bind them?, Adv.Theor.Math.Phys. 

9 (2005) 667-701, hep-th/0408106 . 

[37] H. Elvang, R. Emparan, D. Mateos, and H. S. Reall, 
Supersymmetric black rings and three-charge supertubes, 
Phys.Rev. D71 (2005) 024033, hep-th/0408120 . 

[38] J. P. Gauntlett and J. B. Gutowski, General concentric 
black rings, Phys.Rev. D71 (2005) 045002, 
hep-th/0408122. 

[39] K. Clough, P. Figueras, H. Finkel, M. Kunesch, E. Lim, 
and S. Tunyasuvunakool, to appear. 

[40] L. Lehner and F. Pretorius, Black Strings, Low 
Viscosity Fluids, and Violation of Cosmic Censorship, 
Phys.Rev.Lett. 105 (2010) 101102, arXiv: 1006.5960 . 



